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1.  INTRODUCTION 

Geographical,  seasonal,  and  day-to-day  variations  in  the  vertical  distributions 
of  atmospheric  density  and  temperature  must  be  considered  in  the  design  and 
operation  of  reentry  vehicles.  The  atmospheric  drag  coefficients  for  a  given  vehicle 
are  a  function  of  the  Mach  number,  which  is  directly  related  to  ambient  tempera¬ 
ture.  Variations  in  the  distribution  of  atmospheric  density  affect  the  predicted 
deceleration  and  range  of  free-fall  bombs  and  ballistic  missiles  that  have  high  for¬ 
ward  velocities.  Thus,  the  purpose  of  this  report  is  to  develop  statistical  presenta¬ 
tions  that  can  be  used  to  determine  the  influence  of  density  and  temperature  dis¬ 
tributions  on  the  trajectories  of  reentry  vehicles  at  altitudes  between  the  surface 
and  60  km. 

This  work  is  pari  of  a  continuing  effort  to  compile,  analyze,  and  present  in¬ 
formation  on  the  distribution  of  thermodynamic  properties  of  the  atmosphere  in  a 
form  suitable  for  use  in  the  design  and  operation  of  aerospace  vehicles.  An  econo¬ 
mical  method  of  estimating  the  effect  of  density  and. temperature  on  reentry  vehicles 
by  the  use  of  interlevel  coefficients  of  correlation  of  density  with  density  and  temp¬ 
erature  with  temperature  is  described.  A  technique  is  also  presented  for  obtaining 
estimates  of  these  interlevel  relationships  in  regions  of  sparse  data. 


(Received  for  publication  13  May  1980). 
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In  evaluating  the  accuracy  of  a  proposed  missile  system,  estimates  must  be 
made  of  the  effect  of  day-to-day  variations  in  density  on  the  impact  point.  The 
most  direct  approach  to  this  problem  is  to  calculate  individual  points  of  impact  in 
each  target  area  using  a  representative  sample  of  observed  density  profiles  from 
each  region.  This  task  could  be  time  consuming,  especially  if  several  types  of 
aerodynamic  designs  and  various  reentry  angles  are  being  investigated.  A  simpler 
approach  is  possible  if  the  effect  of  density  at  each  level  is  known  to  be  independent 
or  nearly  independent  of  the  effects  at  other  levels,  as  it  is  believed  to  be  in  most 
cases.  The  required  representative  sample  of  soundings  can  then  be  summarized 
statistically  and  presented  in  an  array  so  that  only  a  simple  arithmetical  calcula¬ 
tion  is  required  to  find  the  missile  range  and  its  standard  deviation. 

Arrays  of  means  and  standard  deviations  of  density  and  temperature  at  2 -km 
intervals  of  altitude  from  the  surface  to  60  Km,  together  with  interlevel  coefficients 
of  correlation  of  density  with  density  and  temperature  with  temperature,  are  pre¬ 
sented  for  the  mid-season  months  at  four  locations.  The  same  type  of  data  are  also 
presented  for  altitudes  between  26  and  60  km  at  six  additional  locations.  Both  sets 
of  data,  ranging  in  latitude  from  8°S  to  64°  N,  are  in  a  format  that  can  be  readily 
used  to  determine  the  integrated  effect  of  density  and  temperature  variations  on  the 
trajectories  of  reentry  vehicles. 

The  integrated  effect  (E)  of  mean  monthly  density  on  the  trajectory  and  impact 
point  of  a  missile  can  be  determined  for  a  specific  location  by  computer-simulated 
"flights"  through  mean  monthly  or  seasonal  density  profiles,  if  the  proper  influence 
coefficients  (c.)  for  the  missile  at  various  levels  are  given: 

E  *  I  c.p.  .  (1) 


where  p.  is  the  mean  monthly  density  at  the  i'th  level.  The  integrated  standard 
deviation  of  the  range,  or  deceleration  due  to  day-to-day  variations  from  the 
mean  seasonal  or  monthly  density  profiles,  can  be  obtained: 


I..  c.  a.  R..  c.  a.  , 
«J  *  >  JJ  J  J 


t’2) 


2 

where  a.nt  is  the  integrated  variance  for  all  layers  being  considered,  c^  and  c.  are 
influence  coefficients  at  levels  i  and  j,  a.  and  a.  are  the  standard  deviations  of 
density  at  levels  i  and  j,  and  R. ^  is  the  coefficient  of  correlation  between  densities 
at  levels  i  and  j. 
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Previous  statistical  presentations  of  the  vertical  density  distributions  suitable 

for  reentry  studies  were  prepared  for  levels  up  to  30  km  from  radiosonde  observa- 
1  2 

tions  by  Cole  and  Court  and  by  Cole  and  Nee.  Investigations  at  altitudes  above 
30  km  had  been  hampered  by  the  lack  of  adequate  data.  Only  during  the  past  few 
years  have  sufficient  observations  become  available  from  Meteorological  Rocket 
Network  (MRN)  stations  for  such  investigations. 


2.  DATA  SOURCES  AND  ACCURACY 
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The  statistical  arrays  of  the  means,  standard  deviations,  and  interlevel  coef¬ 
ficients  of  correlation  of  density  and  temperature  for  altitudes  up  to  60  km  are  based 
on  radiosonde  and  MRN  observations  taken  at  the  10  locations  shown  in  Table  1.  At 
four  of  these  sites,  radiosonde  and  rocketsonde  observations  taken  within  a  few 
hours  of  each  other  were  combined  to  provide  individual  density  and  temperature 
profiles  from  the  surface  to  60  km.  Rocketsondes  without  accompanying  radio¬ 
sondes  were  used  at  the  six  remaining  sites  to  develop  statistical  properties  of  the 
density  and  temperature  distributions  between  26  and  60  km.  Only  observations 
that  were  taken  at  least  3  days  apart  were  used  in  an  attempt  to  reduce  the  effects 
of  persistence  on  the  statistical  properties  of  the  elements  being  studied.  It  has 

3 

been  shown  by  Durst  that  the  coefficients  of  correlation  (R)  for  meteorological 
elements  observed  at  various  time  intervals  at  levels  between  5  and  20  km  follow 
the  law  R  =  e  ai",  where  a  equals  6.  9  X  10  6  sec  1  and  T  is  measured  in  seconds. 

In  a  3-day  interval,  the  coefficients  of  correlation  between  observations  are  gener¬ 
ally  less  than  0.  2,  approaching  independence. 

Temperatures  were  obtained  at  2 -km  intervals  of  geopotential  altitude  below 
30  km  from  individual  radiosonde  observations  by  linear  interpolation  between  values 
at  standard  and  significant  pressure  levels  whose  heights  had  been  computed  hypso- 
metrically.  Densities  were  interpolated  logarithmically.  Errors  introduced  by 

logarithmic  interpolation  of  density  are  unimportant,  generally  <  1  percent  at  radio- 
2 

sonde  heights.  The  root -mean -square  (rms)  observational  errors  in  radiosonde 
temperatures  vary  linearly  with  altitude  from  1°K  at  the  surface  to  from  2  to  5°K 
at  30  km,  and  the  estimated  rms  errors  in  densities  derived  from  these  measure- 

4 

ments  range  from  0.2  percent  at  the  surface  to  1.0  percent  at  30  km. 

1.  Cole,  A.  E. ,  and  Court,  A.  (1962)  Density  Distribution,  Interlevel  Correlations, 

and  Variation  with  Winds,  AFCRL-TR-82 -815.  1 

2.  Cole,  A.  E.,  and  Nee,  P.  F.  (1965)  Correlations  of  Temperature,  Pressure, 

and  Density  to  30  km,  AFCRL-TR-65-43,  AD  612651. 

3.  Durst,  B.  A.  (1954)  Variation  of  Wind  with  Time  and  Distance,  Geophysical 

Memoirs  No.  93,  British  Meteorological  Gffice. 

4.,  Meteorological  Group,  Range  Commanders  Council  (1977)  Meteorological 

Data  Error  Estimates,  Document  110-77,  White  Sands  Missile  &ange,  NM. 
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Temperatures  and  densities  were  obtained  from  MRN  observations  at  geo¬ 
metric  altitudes  between  26  and  60  km.  The  estimated  rms  observational  errors 
in  rocketsonde  temperatures  vary  linearly  with  altitude  from  1.  5°K  at  25  km  to 
3.5°K  at  70  km,  and  the  rms  errors  in  densities  derived  from  rocketsonde  measure- 

4 

ments  are  estimated  to  vary  from  1  percent  at  30  km  to  5  percent  at  60  km. 


Table  1.  Observational  Sites 


Station 

Location 

Altitudes 

Period  of  Record 

Ascension  Island 

8°S, 

14°W 

Surface  to  60  km 

1969-1976 

Kwajalein 

94N, 

Surface  to  60  km 

1969-1976 

Wallops  Island 

38°N. 

75°  W 

Surface  to  60  km 

1969-1976 

Churchill 

59°N, 

94°  W 

Surface  to  60  km 

1969-1976 

Fort  Sherman 

9°N, 

80°W 

26  to  60  km 

1969-1976 

Barking  Sands 

22°N, 

160°W 

26  to  60  km 

1969-1976 

Cape  Kennedy 

28°N, 

80°W 

26  to  60  km 

1969-1976 

White  Sands 

32“N, 

106°W 

26  to  60  km 

1969-1976 

Primrose  Lake 

55°N, 

110°W 

26  to  60  km 

1969-1976 

Poker  Flats 

64<1N, 

146°W 

26  to  60  km 

1969-1976 

The  observed  rms  variations  (Oq)  around  the  monthly  means  contain  the  true 
rms  variability  due  to  changes  in  synoptic  conditions  (o^,)  and  the  rms  observational 
error  (cte ).  If  the  true  variability  and  observational  errors  are  assumed  to  be 
independent,  the  observed  rms  variability  is  given  by: 

“o  ■  V  4  *  4  •  <3> 

Consequently,  the  effect  of  observational  errors  must  be  carefully  evaluated  to 
determine  how  much  of  the  variability  that  is  indicated  by  the  uncorrected  soundings 
is  due  to  synoptic  changes  in  weather  patterns.  These  errors  have  less  effect  on 
the  mean  monthly  values  of  temperature  and  density  given  in  the  appendices,  be¬ 
cause  the  rms  error  of  the  mean  monthly  values  is  equal  to  the  rms  value  of  the 
error  of  an  observation  divided  by  the  square  root  of  the  number  of  independent 
observations  used  for  computing  the  monthly  means. 
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3.  TEMPERATURE 


Statistical  arrays  of  mean  monthly  temperatures,  standard  deviations,  and  in¬ 
terlevel  coefficients  of  correlation  are  given  in  Appendix  A  for  altitudes  up  to  6 Okm. 
An  indication  of  the  variations  in  the  vertical  distributions  of  mean  monthly  tempera¬ 
tures  with  geographical  location  can  be  obtained  from  the  temperature -altitude 
profiles  for  each  of  the  four  mid-season  months  shown  in  Figure  1  at  Ascension 
Island  (8S),  Wallops  Island  (38N),  and  Churchill  (59N).  The  most  striking  feature 
of  these  comparisons  is  seen  at  the  two  levels  where  the  temperature  profiles  for 
the  three  locations  converge  or  cross.  The  first  such  level  occurs  near  12  km; 
the  second,  between  22  and  26  km.  Below  12  km,  the  temperature  profiles  tend  to 
parallel  each  other,  indicating  that  geographical  differences  in  the  mean  monthly 
temperature -lapse  rates  in  the  troposphere  are  relatively  small.  Locations  that 
have  the  warmest  temperatures  below  12  km  have  the  coldest  temperatures  between 
12  and  24  km.  The  largest  geographical  temperature  variations  occur  at  the  sur¬ 
face  in  January,  April,  and  October,  and  at  16  km  in  July.  Above  24  km,  the 
variability  of  the  profiles  with  latitude  is  largest  in  October  and  January. 

Seasonal  changes  in  the  temperature -altitude  profile  at  each  of  the  three  loca¬ 
tions  are  illustrated  in  Figure  2.  Differences  in  mean  monthly  temperatures  are 
largest  at  the  surface  and  smallest  near  12  km.  The  vertical  temperature  gradients 
in  the  troposphere  are  negative  and  remain  relatively  constant  throughout  the  year. 
The  mean  monthly  vertical  gradients  in  the  stratosphere  are  positive  in  all  months 
at  Ascension  Island  and  Wallops  Island;  at  Churchill,  they  are  slightly  negative  or 
isothermal  between  10  and  28  km  in  January,  October,  and  April  and  slightly 
positive  in  July. 

Standard  deviations  of  the  day-to-day  variations  around  the  mean  monthly 
temperatures  at  levels  between  the  surface  and  60  km  are  tabulated  in  Appendix  A 
and  are  plotted  in  Figure  3  for  January  and  July  at  Ascension  Island,  Wallops  Island, 
and  Churchill,  These  variations,  due  to  day-to-day  changes  in  weather  patterns, 
are  largest  in  January  and  at  high  latitudes,  whereas  in  July  the  variations  are 
approximately  the  same  at  all  latitudes.  The  magnitude  of  the  seasonal  difference 
in  the  standard  deviations  of  temperature  at  these  altitudes  increases  with  latitude. 
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Figure  1.  Latitudinal  Changes  in  the  Temperature-Altitude  Profiles  for  the  Mid- 
Season  Months  at  Ascension  Island,  Wallops  Island,  and  Churchill 
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Figure  3.  Day-to-Day  Variability 
Around  Mean  Monthly  Temperatures 
for  January  and  July  at  Ascension 
Island,  Wallops  Island,  and 
Churchill 


ST4JCAR0  DEVIATION  (PERCENT  OF  MEAN)  Of  DENSITY 


The  day-to-day  variations  around  mean  monthly  densities  are  plotted  vs 
altitude  in  Figure  6  as  coefficients  of  variation  (100  X  SD/mean)  for  the  mid-season 
months  at  Ascension  Island,  Wallops  Island,  and  Churchill.  These  variations 
increase  markedly  with  latitude,  particularly  during  winter.  In  July,  however, 
they  are  essentially  the  same  at  all  latitudes. 

4.  DENSITY 

Arrays  of  mean  monthly  densities,  standard  deviations,  and  interlevel  correla¬ 
tions  for  auitudes  up  to  60  km  are  given  in  Appendix  B.  The  mean  densities  for  each 
of  the  mid-season  months  at  Ascension  Island,  Wallops  Island,  and  Churchill  are 
plotted  versus  altitude  in  Figure  4  as  percentage  departures  from  the  U. S.  Standard 
Atmosphere,  1976.  The  individual  profiles  cross  or  converge  near  8  km  and  be¬ 


tween  22  and  26  km.  Both  are  levels  of  minimum  density  variability.  There  is  an 
isopycnic  level  near  8  km  at  which  mean  monthly  densities  depart  from  standard  by 
no  more  than  1  or  2  percent,  regardless  of  location  and  season.  However,  between 
22  and  26  km  there  is  a  marked  seasonal  variability,  even  though  there  is  very  little 
latitudinal  variability. 

Seasonal  changes  n  the  mean  monthly  density  profiles  at  each  of  the  three  loca¬ 
tions  are  shown  in  Figure  5.  The  minimum  seasonal  variability,  1  or  2  percent, 
occurs  at  8  km  and  the  maximum  occurs  above  60  km. 


5.  COESA  (1976)  U.  S.  Standard  Atmosphere,  1976,  GFO,  Washington,  D.  C. 
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Figure  4.  Lutltudinal  Changes  in  tho  Density-Altitude  Profiles  for  the  Mid-Season 
Months  ut  Ascension  Island,  Wallops  Island,  and  Churchill 
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Figure  5.  Seasonal  Differences  in  the  Density -Altitude  Profiles  at  Ascension  Island, 
Wallops  Island,  and  Churchill 
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The  day-to-day  variations  around  mean  monthly  densities  are  plotted  vs 
altitude  in  Figure  6  as  coefficients  of  variation  (lOOxSD/mean)  for  the  mid-season 
months  at  Ascension  Island,  Wallops  Island,  and  Churchill.  These  variations 
increase  markedly  with  latitude,  particularly  during  winter.  In  July,  however, 
they  are  essentially  the  same  at  all  latitudes. 
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5.  INTERLEVEL  CORRELATION 

The  manner  in  which  the  correlation  between  densities  at  two  levels  decreases 
(or  decays)  with  increasing  separation  between  the  levels  is  an  example  of  the 
general  problem  of  correlation  decay.  Similar  correlation  decay  is  found  for  most 
meteorological  elements  as  the  horizontal  or  vertical  distance  between  the  points  of 
observation  or  the  time  interval  between  observations  increases.  As  yet,  no  fully 
satisfactory  description  of  the  decay  rate,  based  on  fundamental  properties  or 
assumptions,  is  available.  Consequently,  many  empirical  models  that  are  valid 
for  specific  elements  over  restrictive  ranges  have  been  proposed. 

Profiles  of  coefficients  of  correlation,  R,  of  surface  density  with  density  at 
other  altitudes  are  shown  in  Figure  7  for  each  of  the  mid-season  months  at  Kwajalein, 
Wallops  Island,  and  Churchill.  During  all  four  seasons  and  at  all  three  locations  the 
correlation  between  surface  density  and  density  at  other  altitudes  decreases  rapidly 
as  the  vertical  distance  between  levels  increases,  approaching  zero  by  8  or  10  km. 
The  correlation  then  generally  remains  near  zero  or  slightly  negative  at  altitudes 
between  10  and  60  km.  The  low  correlation  between  the  surface  density  and  the 
density  at  levels  between  10  and  60  km  indicates  that  very  little  information  on 
deviations  from  the  mean  monthly  density  profiles  between  10  and  60  km  can  be  ob¬ 
tained  from  observations  of  surface  density. 

Profiles  of  coefficients  of  correlation  of  the  density  at  26  km  with  density  at 
higher  altitudes  up  to  60  km  are  shown  for  each  of  the  mid-season  months  for  loca¬ 
tions  in  the  tropics  (Figure  8),  middle  latitudes  (Figure  9),  and  high  latitudes 
(Figure  10).  These  figures  indicate  that  at  middle  and  high  latitudes  the  interlevel 
density  correlations  above  26  km  decay  more  rapidly  in  winter  than  in  the  other 
seasons,  and  that  seasonal  variations  in  the  decay  rates  in  the  tropics  are  relatively 
small.  In  April,  July,  and  October  the  26-km  densities  at  most  locations  are  posi¬ 
tively  correlated  with  the  densities  at  nearly  all  levels  between  25  and  60  km.  At 
several  stations,  data  on  which  to  base  these  profiles  of  interlevel  correlations  were 
not  available  for  all  four  mid- season  months.  As  a  result  the  April  profile  for 
Ascension  Island  is  not  shown  in  Figure  8,  and  the  Primrose  Lake  and  Poker  Flats 
profiles  for  July  and  October,  respectively,  are  not  shown  in  Figure  10. 

Profiles  of  coefficients  of  correlation,  R,  of  surface  temperatures  with  temper¬ 
atures  at  other  altitudes  are  shown  in  Figure  11  for  each  of  the  mid-season  months 
at  Ascension  Island,  Kwajalein,  Wallops  Island,  and  Churchill.  At  most  locations, 
the  correlation  between  surface  temperatures  and  temperatures  at  other  altitudes 
decreases  rapidly  with  increasing  altitude,  reaching  a  minimum  or  becoming  nega¬ 
tive  between  12  and  16  km  and  then  remaining  near  zero,  plus  or  minus  0.  3,  from 
20  to  60  km.  There  is  a  pronounced  negative  correlation  between  the  surface 
temperature  at  Wallops  Island  and  temperatures  at  altitudes  between  10  and  20  km. 
Values  of  -0.4  to  -0.6  occur  at  these  levels  in  January,  April,  and  October. 
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Figure  8.  Vertical  Profiles  of  Interlevel  Coefficients  of  Correlation  of  Density 
at  26  km  With  Density  at  Higher  Altitudes  up  to  60  km  for  the  Mid-Season 
Months  at  Ascension  Island,  Fort  Sherman,  and  Kwajalein 
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Figure  9.  Vertical  Profiles  of  Interlevel  Coefficients  of  Correlation  of  Density  at 
26  km  With  Density  at  Higher  Altitudes  up  to  60  km  for  the  Mid-Season  Months  at 
Barking  Sands,  Cape  Kennedy,  White  Sands,  and  Wallops  Island 
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Figure  10.  Vertical  Profiles  of  Irrterlevel  Coefficients 
26  km  With  Density  at  Higher  Altitudes  up  to  60  km  for 
Primrose  Lake,  Churchill,  and  Poker  Flats 
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Figure  11.  Vertical  Profiles  of  Interlevel  Coefficients  of  Correlation  of  Surface 
Temperature  With  Temperature  at  Other  Altitudes  up  to  60  km  for  the  Mid-Season 
Months  at  Ascension  Island.  Kwajalein,  Wallops  Island,  and  Churchill 
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Individual  profiles  of  the  correlations  between  temperatures  at  26  km  and  those 
at  other  altitudes  between  26  and  60  km  are  plotted  in  Figures  12,  13,  and  14  for 
locations  at  low,  middle,  and  high  latitudes,  respectively.  The  July  profile  for 
Primrose  Lake  is  not  shown  in  Figure  14  because  only  seven  observations  were 
available  for  the  development  of  the  interlevel  correlations.  The  rate  of  decay  of 
the  interlevel  correlations  with  increasing  vertical  separation  between  levels  is 
approximately  the  same  at  all  latitudes  in  July,  and  positive  correlation  exists  be¬ 
tween  the  values  at  26  km  and  those  of  other  altitudes  up  to  60  km.  In  January  the 
rates  of  decay  at  middle  and  high  latitudes  are  considerably  greater  than  those  in 
the  tropics,  and,  except  in  the  tropics,  negative  correlations  occur  at  most  loca¬ 
tions  when  levels  are  separated  by  more  than  13  to  20  km. 

The  individual  arrays  of  interlevel  correlations  of  temperatures  and  densities 
from  26  to  60  km,  given  in  Appendices  A  and  B  for  locations  in  tropical  areas 
(Ascension  Island,  Fort  Sherman,  and  Kwajalein),  were  combined  to  produce  one 
set  each  of  temperature  and  density  matrices,  Tables  2  and  3  respectively,  for 
tropical  areas.  This  same  procedure  was  followed  for  middle  latitudes  (Barking 
Sands,  Cape  Kennedy,  White  Sands,  and  Wallops  Island),  Tables  4  and  5,  and  for 
high  latitudes  (Primrose  Lake,  Churchill,  and  Poker  Flats),  Tables  6  and  7.  Com¬ 
bining  the  individual  arrays  increased  the  sample  size  and  smoothed  out  many  of  the 
irregularities  found  in  the  profiles  of  R  that  are  ba§ed  on  values  from  one  station. 

The  combined  or  average  11-values  of  density  are  plotted  at  2-km  intervals  between 
altitudes  of  26  and  60  km,  using  semi-log  coordinates  (Figures  15,  16,  and  17)  for 
low,  middle,  and  high  latitudes,  respectively.  11-values  cf  temperature  are  pre¬ 
sented  in  the  same  manner  in  Figures  18,  19,  and  20.  The  26-km  level  was  selected 
as  the  lower  limit  because  it  is  near  the  top  of  observed  radiosonde  temperature  and 
density  profiles  that  are  taken  on  a  routine  basis  in  most  areas  of  the  world.  Conse¬ 
quently,  temperatures  and  densities  observed  by  radiosondes  at  26  km  can  be  used 
with  the  interlevel  correlations  of  temperatures  and  densities  shown  in  Tables  2 
through  7  to  extrapolate  these  elements  to  higher  altitudes. 
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Figure  14.  Vertical  Profiles  of  Interlevel 
Coefficients  of  Correlation  of  Temperature 
at  26  km  With  Temperature  at  Higher 
Altitudes  up  to  60  km  for  the  Mid-Season 
Months  at  Primrose  Lake,  Churchill,  arid 
Poker  Flats 


Table  2a.  Average  January  Interlevel  Corre¬ 
lations  of  Temperature  for  the  Tropics 
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Table  2c.  Average  July  Interlevel  Corre¬ 
lations  of  Temperature  for  the  Tropics 
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Table  2b.  Average  April  Interlevel  Corre¬ 
lations  of  Temperature  for  the  Tropics 
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Table  2d.  Average  October  Interlevel  Corre¬ 
lations  of  Temperature  for  the  Tropics 
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Multiply  tabular  values  by  0,01  to  obtain  correlation  coefficients. 
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Table  3d.  Average  October  Interlevel  Corre 
lations  of  Density  for  the  Tropics 
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Multiply  tabular  values  by  0.01  to  obtain  correlation  coefficients. 
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Table  4c.  Average  July  Interlevel  Corre¬ 
lations  of  Temperature  for  Middle  Latitudes 
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Multiply  tabular  values  by  0.  01  to  obtain  correlation  coefficients 


Table  5a.  Average  January  Interlevel  Corre¬ 
lations  of  Density  for  Middle  Latitudes 
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Table  5c.  Average  July  Interlevel  Corre¬ 
lations  of  Density  for  Middle  Latitudes 
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Multiply  tabular  values  by  0.  01  to  obtain  correlation  coefficients, 
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Table  6c.  Average  July  Interlevel  Corre- 
lations  of  Temperature  for  High  Latitudes 
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58  57  55  47  7  1  57  50  63  57  3  3  3  4  3  2  6  5  7  2  80  8  1  86 

60  4  6  47  32  61  49  41  5*  56  30  36  3  6  65  67  73  7  0  80  >! 


Table  6d.  Average  October  Interlevel  Corre 
lations  of  Temperature  for  High  Latitudes 


_  76  26  30  32  34  36  38  40  42  44  46  48  50  52  54  56  58 

28  85*  ~ 

30  82  89 

32  75  80  93 

34  62  73  85  92 

36  54  64  76  85  8« 

38  51  60  71  80  85  87 

40  49  53  64  74  77  79  93 

42  31  36  44  56  66  69  82  86 

44  24  32  37  51  62  65  77  80  91 

46  16  25  23  3*  48  52  63  64  82  *8 

48  13  23  22  34  44  49  57  58  73  81  92 

08  55  15  24  33  32  47  51  67  7l  7<»  87 

07  09  09  |7  29  21  41  48  63  65  71  77  91 

54  03  08  08  20  2?  23  33  3«  48  48  55  59  69  82 

56  02  -01  -01  09  14  12  22  31  37  35  39  47  53  67  82 

58  -15  -11  -12  -03  03  09  14  16  21  28  32  41  46  56  71  89 

09  -28  -30  -13  -06  -30  -26  15  17  28  33  37  30  50  62  85 


Multiply  tabular  values  by  0.01  to  obtain  correlation  coefficients 


Table  7a.  Average  January  Interlevel  Corre¬ 
lations  of  Density  for  High  Latitudes 


km 

2C 

28 

30 

32 

34 

36 

38 

40 

42 

44 

46 

48 

50 

52 

54 

56 

58 

28 

55' 

30 

52 

97 

32 

8C 

94 

97 

34 

71 

89 

93 

98 

3* 

72 

84 

18 

94 

98 

38 

Cl 

74 

80 

86 

92 

96 

40 

45 

C3 

72 

79 

86 

•»2 

97 

42 

41 

58 

65 

72 

80 

87 

53 

*7 

44 

37 

52 

60 

CC 

74 

81 

90 

93 

58 

46 

30 

4$ 

52 

58 

CC 

73 

82 

87 

93 

97 

48 

22 

37 

43 

50 

57 

65 

75 

cl 

8« 

A4 

98 

50 

14 

28 

85 

40 

47 

55 

CC 

73 

82 

88 

94 

98 

52 

10 

24 

24 

32 

26 

44 

53 

60 

71 

78 

86 

91 

95 

54 

01 

IS 

2  i 

23 

29 

36 

51 

53 

64 

73 

81 

-»o 

96 

95 

SC 

03 

12 

20 

19 

23 

28 

37 

46 

58 

67 

78 

87 

93 

*»3 

99 

58 

•04 

o; 

15 

15 

20 

26 

*6 

46 

58 

65 

76 

84 

90 

91 

96 

99 

80 

•01 

12 

19 

1“ 

24 

27 

33 

45 

57 

62 

6* 

77 

83 

35 

91 

94 

98 

Table  7b.  Average  April  Interlevel  Corre¬ 
lations  of  Density  for  High  Latitudes 


km 

26 

28 

30 

32 

34 

36 

38 

40 

42 

44 

46 

48 

50 

52 

~54 

56 

58 

28 

92 

30 

8) 

93 

32 

76 

87 

94 

34 

69 

80 

85 

94 

36 

59 

69 

76 

84 

93 

38 

S3 

60 

65 

73 

84 

93 

40 

45 

52 

54 

64 

76 

88 

95 

42 

39 

41 

42 

48 

63 

75 

86 

93 

44 

32 

34 

34 

37 

50 

64 

77 

84 

46 

31 

33 

32 

33 

44 

58 

70 

78 

88 

95 

48 

28 

27 

25 

25 

35 

49 

63 

72 

84 

92 

97 

50 

24 

21 

19 

«' 

27 

41 

54 

64 

79 

89 

94 

98 

52 

18 

15 

13 

10 

19 

33 

47 

57 

74 

86 

93 

96 

98 

54 

19 

16 

14 

09 

18 

31 

44 

52 

70 

83 

91 

91 

97 

99 

56 

17 

14 

11 

07 

14 

26 

38 

46 

64 

78 

87 

93 

95 

98 

99 

58 

12 

0b 

04 

•  02 

06 

20 

33 

42 

62 

78 

85 

90 

>4 

97 

97 

98 

60 

00 

05 

03 

02 

11 

26 

39 

49 

65 

76 

86 

90 

93 

93 

96 

97 

99 

Table  7c.  Average  July  Interlevel  Corre¬ 
lations  of  Density  for  High  Latitudes 


km 

26 

28 

30 

32 

34 

36 

36 

40 

42 

44 

46 

48 

50 

>2 

54 

55 

58 

28 

94* 

30 

88 

90 

32 

81 

88 

86 

34 

71 

80 

82 

94 

36 

51 

61 

66 

88 

80 

38 

46 

59 

52 

76 

72 

83 

40 

42 

54 

44 

68 

68 

74 

88 

42 

26 

42 

38 

64 

64 

68 

78 

8* 

44 

13 

28 

26 

54 

55 

70 

74 

80 

>2 

46 

n 

25 

22 

52 

52 

72 

78 

77 

80 

*»2 

48 

07 

20 

20 

SO 

51 

73 

72 

74 

76 

85 

®4 

50 

05 

20 

25 

4* 

45 

55 

68 

66 

70 

78 

52 

92 

52 

06 

15 

19 

44 

42 

62 

68 

68 

75 

80 

30 

88 

96 

54 

08 

18 

18 

48 

44 

67 

74 

72 

7* 

84 

85 

88 

91 

96 

56 

00 

10 

10 

42 

44 

65 

72 

75 

82 

88 

*0 

86 

85 

87 

96 

58 

08 

12 

10 

40 

45 

60 

62 

72 

82 

*1 

88 

S3 

76 

SO 

89 

96 

60 

06 

IS 

10 

16 

42 

62 

62 

•  1 

83 

92 

90 

83 

73 

'i 

87 

96 

99 

Table  7d.  Average  October  Interlevel  Corre¬ 
lations  of  Density  for  High  Latitudes 


km 

26 

28 

30 

32 

34 

36 

38 

40 

= 

42 

44 

46 

48 

50 

52 

54 

56 

58 

28 

84 

30 

73 

90 

32 

64 

80 

*2 

34 

60 

73 

86 

94 

36 

51 

65 

76 

87 

93 

38 

44 

55 

6? 

83 

91 

96 

40 

38 

48 

63 

7* 

85 

?2 

42 

36 

44 

59 

75 

52 

91 

95 

98 

44 

36 

42 

55 

71 

78 

89 

93 

95 

56 

46 

3! 

38 

50 

66 

72 

90 

93 

96 

93 

48 

31 

3? 

50 

64 

70 

81 

*7 

9J 

94 

97 

98 

50 

26 

33 

45 

S» 

66 

78 

84 

88 

92 

95 

97 

99 

52 

27 

34 

45 

5? 

65 

76 

82 

56 

90 

94 

96 

98 

>9 

54 

25 

33 

44 

56 

Cl 

74 

79 

83 

87 

91 

94 

97 

98 

9* 

56 

25 

32 

43 

55 

58 

■s 

76 

<1 

65 

69 

92 

95 

97 

98 

99 

58 

25 

34 

43 

56 

59 

71 

77 

83 

87 

50 

93 

96 

97 

98 

98 

99 

60 

24 

35 

43 

54 

(0 

70 

77 

82 

65 

89 

92 

95 

06 

96 

97 

98 

99 

Multiply  tabular  values  by  0.01  to  obtain  correlation  coefficients. 


The  decay  of  the  interlevel  correlations  of  both  temperatures  and  densities 
between  values  at  26  km  and  those  at  higher  altitudes.  Figures  15  through  20,  can 
be  defined  by  a  simple  expression  of  the  type: 

R  =  exp  (-ah)  , 

where  a  is  a  constant  and  h  is  the  vertical  separation  in  kilometers.  Howev  r, 
several  of  the  curves  were  better  defined  by  modified  versions  of  the  same  ex¬ 
pression: 

2 

R  =  exp  (-ah  )  (5 

and 

R  =  exp  (-a-^h)  .  (6 

It  was  found  that  the  value  of  the  constant  "a"  in  these  equations,  which  are  gener¬ 
ally  valid  for  thicknesses  of  12  to  20  km,  changed  with  season  and  latitude,  as 
shown  in  Table  8  and  in  Figures  15  through  20. 


Table  8.  Seasonal  and  Latitudinal  Changes  in  "a"  (used  for  defining  cor¬ 
relation  decay) 


Density 

High  Latitude 

Middle  Latitude 

Low  Latitude 

Jan 

0.  0036h2 

0.  022h2 

0. 195h 

Apr 

0.  06  Oh 

0.  115h 

0.  14  Oh 

July 

0.  006h2 

0.  055h 

0.  130h 

Oct 

0.  06  Oh 

0.  IlOh 

0. 140h 

Temperature 

Jan 

9 

0.  01  Oh 

0.  020h2 

0.  370Vh" 

Apr 

0. 020h2 

0.  130h 

0.  350-^h" 

July 

0.  036h 

0.  07  Oh 

0.  270  Vh~ 

Oct 

0.  0045h2 

0.  20^ 

0.  220  ^h~ 
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If  density  observations  on  a  specific  day  extend  only  to  an  altitude  of  26  km,  the 
usual  peak  altitude  of  a  radiosonde,  ar  estimated  density  profile  can  be  obtained 
from  28  to  60  km  by  using  the  appropriate  monthly  means,  standard  deviations,  and 
interlevel  correlations  in  the  general  expression: 

K=  h  +  R^-(pi  (7> 

A  1 

where  p„  is  the  estimated  density  at  the  desired  level,  p.  is  the  density  at  the  lower 
^  * 
level,  p2  is  the  density  at  the  upper  level,  Oj  is  the  standard  deviation  of  Pj.  Og 

is  the  standard  deviation  of  p2»  and  R  is  the  coefficient  of  correlation  of  density 

between  the  two  levels. 

The  effect  of  density  on  the  trajectory  and  impact  point  of  a  missile  can  be 
estimated  through  use  of  Eqs.  (1)  and  (2).  The  mean  effect,  E,  can  be  determined 
from  Eq.  (1)  for  a  particular  location  or  region,  using  the  mean  monthly  densities 
from  the  appropriate  array  in  Appendix  B  and  the  influence  coefficients  for  a  given 
missile  at  the  various  levels  of  interest.  The  standard  deviation  around  the  mean 
monthly  effects  of  density  due  to  day-to-day  variations  in  the  density  profile  can  be 
estimated  from  Eq.  (2).  An  arithmetic  example,  given  in  Appendix  C,  illustrates 
the  use  of  Eqs.  (1)  and  (2)  along  with  the  statistical  arrays  from  Appendix  B. 

6.  VERTICAL  DENSITY  GRADIENTS 

The  largest  changes  of  density  with  altitude  occur  in  layers  that  have  a  low  base 
temperature  and  a  large  positive  vertical  temperature  gradient.  Consequently,  the 
largest  mean  monthly  vertical  density  gradients  below  3C  km  occur  in  the  layer 
immediately  above  the  tropical  tropopause  where  temperatures  increase  rapidly 
withaltitudefromtropospheric and  stratospheric  minima  near  190°K.  Vehicles 
descending  through  this  region  would  encounter  a  44-  to  51 -percent  increase  in 
density  through  a  2-km  layer  (Tables  B1  and  B2).  The  largest  density  gradients 
above  30  km  within  the  stratosphere  and  lower  mesosphere  generally  occur  between 
30  and  40  km  in  arctic  and  subarctic  regions  is?  January,  when  temperatures  are 
lowest  near  30  km  and  increase  with  altitude  up  to  45  or  50  km.  This  can  be  even 
more  significant  in  winter  during  a  "stratospheric  warming"  in  which  the  tempera¬ 
ture  in  the  lower  stratosphere  decreases  while  the  temperature  in  the  upper 
stratosphere  increases  by  30  or  40<>K. 

6.  Kantor.  A.J..  and  Cole,  A.  E.  (1977)  Monthly  90°  N  Atmospheres  and  Hig^- 
Eatitude  Warm  and  Cold  Winter  Stratosphere/ Mesosphere,  AT  GL- TB-»7- 


Mean  monthly  vertical  density  gradients  for  2-km  increments  of  altitude  may 
be  extracted  from  the  statistical  arrays  given  for  10  locations  in  Appendix  B.  More 
importantly,  the  variance  of  density  gradients  associated  with  the  day-to-day  varia¬ 
tions  in  synoptic  weather  conditions  can  be  estimated  (from  the  appropriate  means, 
standard  deviations,  and  coefficients  of  correlation)  using  the  fundamental  relationship; 

a2  =  (Xj  +  f/2  -  2  HcTj  o2  ,  (8) 

2  2 
where  0^  is  the  mean  monthly  density  variance  at  level  1,  <7 2  is  the  variance  at 

level  2,  R  is  the  coefficient  of  correlation  of  density  between  values  at  levels  1 
*  2 

and  2,  and  a  is  the  estimated  variance  around  the  mean  monthly  gradient  between 
these  levels.  A  sample  computation  for  Churchill  in  January  is  given  in  Appendix  D 
for  a  vertical  density  gradient  that  is  equalled  or  exceeded  2.  5  percent  of  the  time. 
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MULTIPLY  MEAN  BY  INDICATED  NEGATIVE  POMER  OF  10 

MULTIPLY  TABULAR  VALUES  BY  0.(1  TO  OBTAIN  CORRELATION  COEFFICIENTS 


Table  132.  Kwajalein— Correlation  of  January  Density  (kg  m  )  From  Surface  to  60  km 
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•*  MULTIPLY  TABULAR  VALUES  3Y  0. 01.1  TO  09TAIN  CORRELATION  COEFFICIENTS 
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•  MULTIPLY  MEAN  BY  IHDICATEO  NEGATIVE  POWER  OF  10 
"  MULTIPLY  TABULAR  VALUES  IT  0.01  TO  OBTAIN  CORRELATION  COEFFICIENTS 


Tabic  B4,  Churchill  — Correlation  of  January  Density  (kg  m  )  K rom  Surface  to  60  km 
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Table  B8.  White  Sands— Correlation  of  l.  lv  -jensuy  (kg  m  )  From  26  km  to  60  km  (Cont> 
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•  MILT  I  PLY  PEAN  3T  inoicateo  negative  POWER  OF  10 
•*  MULTIPLY  TABULAR  VALUES  BY  0.C1  TO  OBTAIN  CORRELATION  COEFFICIENTS 


KM  KILOMETERS  ABOVE  SEA  LEVEL 

MEAN  AVERAGE  OF  OBSERVED  VALUES 

STOV  STANDARD  DEVIATION  OF  VALUES 
IN  PERCENT  OF  MEAN  TINES  10 

N  NUMBER  OF  VALUES  AT  EACH  ALTITUDE 
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•  MULTIPLY  MEAN  BY  INDICATED  NEGATIVE  POWER  OF  10 
••  MULTIPLY  TABULAR  VALUES  BY  0.01  TO  OBTAIN  CORRELATION  COEFFICIENTS 
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Appendix  C 


Sample  Calculation  for  Estimating  the  Effect  of 
Density  on  a  Reentry  Vehicle 


'( 


The  computations  provided  in  this  appendix  illustrate  how  to  estimate  the  effect 
that  the  mean  monthly  density  and  its  day-to-day  variations  around  the  mean  have 
on  the  trajectory  of  a  reentry  vehicle  passing  through  the  region  between  11  and 
5  km  at  Wallops  Island  in  January. 

For  the  purposes  of  this  simplified  example,  the  following  influence  coefficients 
for  the  vehicle  have  been  assumed,  such  that  a  density  of  1  kg  m  3  at  a  specific 
altitude  will  have  the  indicated  retarding  effect  through  a  2-km  layer  of  the  atmos¬ 
phere: 


Layer  (km) 


-3-1 

Influence  (m(kg  m  )  ) 


11-9 

9-7 

7-5 


800 

400 

200 


J 


if 


ti 
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The  average  distance  short  of  a  vacuum  trajectory  is  the  sum  of  the  effects  due  to 
the  mean  monthly  density  at  each  level,  as  determined  from  Eq.  (1)  in  the  intro¬ 
duction  to  the  main  text: 


Level 

(km) 

Wallops  Density 
(kg/m3) 

Influence 
(m  (kg  m-3)-1) 

Distance 

(m) 

10 

0.411 

X  800 

— 

329 

8 

0.524 

X  400 

= 

210 

5 

0.  658 

X  200 

a 

132 

Average  distance  short  of  trajectory 

= 

671  m 

The  integrated  standard  deviation  of  the  deceleration  due  to  day-to-day  variations 
in  the  density  profile  is  determined  from  Eq.  (2)  using  standard  deviations  of  density 
at  levels  1 0,  8,  and  6  km  and  coefficients  of  correlation  between  these  levels  as 
indicated  in  the  array  for  Wallops  Island,  January,  in  Appendix  B. 


Coefficient 

of  Standard  Standard 

Correlation  Deviation  Influence  Deviation  Influence  Distance 


133 

10 


(1.44X  10-2 

X 

800)2 

(7.86X  10"3 

X 

400)2 

(6.58X  10"3 

X 

200)2 

2  X  0.86 

X 

1.44X  10-2 

X 

800 

X 

2  X  -0.  09 

X 

1. 44  X  10"2 

X 

800 

X 

2  X  0.25 

X 

7.  86 X  10“3 

X 

400 

X 

2  2 

■  .  a  =  206  m  ,  o  =  14.4  m,  and  2o  =  29  m. 


2 


7.86X  10'3 

X 

400 

= 

62 

6.58X  10-3 

X 

200 

= 

-  3 

6. 58  X  10'3 

X 

200 

= 

2 

Total 

= 

206  m2 

Based  on  an  assumption  of  a  normal  distribution  of  atmospheric  density,  the  range 
of  the  missile  will  vary  by  less  than  30  m  from  the  mean  monthly  impact  point 
95  percent  of  the  time.  If  the  aim  point  takes  into  consideration  the  mean  monthly 
density  profile,  the  variations  in  the  impact  point  due  to  day-to-day  variations  in 
the  density  profile  will  be  within  ±  30  m  of  the  target  95  percent  of  the  time. 


Appendix  D 


Sample  Calculation  for  Estimating  an  Extreme 
Vertical  Density  Gradient 


The  following  is  ?.n  illustration  of  how  to  estimate  an  extreme  vertical  density 
gradient  that  may  be  encountered  by  an  aerospace  vehicle. 

The  answer  to  the  question  "What  is  the  vertical  density  gradient  that  will  be 
equalled  or  exceeded  2.  5  percent  of  the  time  between  40  and  42  km  near  Churchill 
in  January?",  may  be  obtained  by  using  the  following  expression  taken  from  Eq.  (8): 


a  =  Vffi  +  a2  "  2  RCTia2  ’ 


(Dl) 


2  2 
where  is  the  mean  monthly  density  variance  at  level  1  (40  km),  a ^  is  the  variance 


at  level  2  (42  km),  R  is  the  coefficient  of  correlation  of  density  between  values  at 
40  and  42  km,  and  <7  is  the  standard  deviation  around  the  mean  monthly  gradient 
between  these  levels.  From  Appendix  B  for  Churchill  in  January, 


cr,„  =  13.  8  percent  of  3.  47  X  10"3  or  4.79  X  10'3  kg  m'3 


40 


cr^2  =  4  percent  of  2.  56  X  10  3  or  3.  69  X  10  3  kg  m  3 


Mean  monthly  diff  =  0.  91  X  10  3  or  9.  10  X  10  ^  kg  m  3 


R  =  0.  98 
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From  Eq.  (Dl), 


